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Abstract. We model the hot and dense strongly interacting mater produced in high 
energy heavy ion collisions using relativistic hydrodynamics. Several different sources 
of real photons produced during these collisions are considered and their relative 
importance is assessed. We include contributions from QCD jets, which are allowed 
to loose and gain energy as they proceed through the hot matter. This is treated 
within the AMY formalism. We obtain photon spectra, R\^, and v\ in agreement 
with measurements performed by the PHENIX collaboration. 



1. Introduction 

The measurement of electromagnetic observables in relativistic nuclear collisions has 
constituted a vigorous experimental effort at different facilities around the world for 
several decades, and it is fair to write that this program has lived up to its promises. 
Indeed, photons (real and virtual) carry information complementary to that in hadronic 
observables, as they suffer essentially no final state interaction. One may show [1] 
that the photon emission rate in a finite-temperature medium is directly related to 
the in-medium retarded self-energy of the photon: a valuable quantity. However, 
the measured signal is comprised of several contributions which have to be under 
simultaneous theoretical control before significant claims can be substantiated. 



2. Photon production 

2.1. Non-thermal processes 

As in nucleon-nucleon interactions, nucleus-nucleus collisions will produce prompt 
photons during the early collisions of the nuclear constituents. The prompt (direct 
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+ fragmentation) photon spectrum can be calculated in pp collisions at next-to- 
leading order [2]. The direct contribution comprises both Compton and annihilation 
components: 

E d(Tp ™ mpt = V / dx a dx b f a (x a ,Q)f b (x b ,Q)K diT E^^^-(p,Q) 
d 3 p d 3 p 

+ J2 f dx a dx b f a (x a ,Q)f b (x b ,Q) [ ^K iv ^ d ° ab ^ cd (Q)D r/c (z,Q F ) (1) 
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The NLO effects are contained in the K-factors; their list of arguments has been 
omitted for brevity. In principle, the prompt photon spectrum depends on the choice 
of renormalization, factorization (Q), and fragmentation (Qf) scales; in the above the 
first has been set equal to the second. Choosing them to be pr/V2, where pr is the 
photon transverse momentum, yields a photon spectrum in proton-proton collisions 
in good agreement with observation [3]. This can be then be used as a baseline for 
calculations of photon production in nucleus-nucleus collisions. In this case, isospin 
effects are included in the parton distribution function, by using the parton species for 
the proton or for the neutron, as appropriate. Shadowing corrections are also taken into 
account. 

2.2. Thermal processes 

In a thermal environment, whether composed of partons or of hadrons in the confined 
phase, the photons produced are calculated simply by performing an integral over the 
space-time history of the strongly interacting system: 
/V 7 r rl 3 R 7 
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The differential photon production rate, R, is evaluated in the local thermal rest frame, 
and at a given time will consist of partonic and/or hadronic contributions. 

Another process which depends on the temperature is that of jet-thermal photon 
production. This happens when a jet crossing the hot medium undergoes a photon- 
producing interaction with a thermal parton. This can further be divided in two classes, 
those which will benefit from collinear enhancements, and those which will not [3]. All 
of the rates introduced so far need to be integrated using with a realistic time-evolution 
model of the relativistic nuclear collision. We will use a two-dimensional relativistic 
hydrodynamical approach - AZHYDRO - which provides a good global description 
of hadronic observables [6] and, importantly, is not tuned to specifically reproduce 
electromagnetic measurements. 

3. Jet energy loss and gain 



The initial jet energy profiles calculated from QCD are evolved in time according to 
AMY [5]. The evolution of the parton energy profile, P(E,t), is determined by solving 
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Figure 1. Spectra of produced photons at RHIC, for two different centrality classes. 
The different elements of the theoretical calculations are shown separately. Data are 
from the PHENIX collaboration [HE!]- 



a set of coupled Fokker-Planck rate equations which have the following generic form: 



dP(E) 
alt 
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dr(E + u,u) 
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The transition rate for processes where partons of energy E loose energy u is 
dY{E,uf) I dudt. The energy gain channels are included through the possibility for uj 
to be negative. Also, AMY smoothly handles the crossover between the Bethe-Heitler 
and Landau- Pomeranchuk-Migdal regimes [7j. A free parameter of the model is a s , the 
strong coupling, which is adjusted to reproduce R A \- 



4. Results and conclusions 



The calculated photon spectra, together with their different components, are shown in 
Figure [TJ The reader may judge of the good agreement between theory and data. No 
parameters have been tuned to generate these figures. Two other useful representations 
of the experimental data are displayed in Figure [2] The first (left panel) is the ubiquitous 
Raa variable, while the second is a plot of the photon azimuthal asymmetry, v\. In both 
panels, the data does not yet permit an accurate determination of model parameters 
and ingredients. However, the fact that R AA < 1 at high values of px is a consequence 
of the fragmentation photons being affected by the energy loss of the fragmenting jets 
(as well as by an isospin effect [3]): should this trend be confirmed experimentally, a 
link would exist between the nuclear modification factor of photons and that of strongly 
interacting particles also borne out of fragmentation. 

Turning to vj, one observes that the azimuthal anisotropy associated with the jet- 
plasma photons is indeed negative, as originally suggested [ID] . However, the net value 
of v J appears to be essentially zero for px > 4 GeV. The reason for this is two-fold. 
First, the geometric anisotropy of the medium in the realistic 2D+1 model is rather 
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Figure 2. The nuclear modification factor of produced photons at RHIC (left 
panel), and the photon azimuthal anisotropy (right panel), for two different centrality 
classes. For R? AA , cases with and without jet-plasma photons are shown. The different 
components of the photon contribution and their azimuthal anisotropy are shown on 
the right. Data are from the PHENIX collaboration. 



small initially, and then proceeds to shrink as time evolves. Second, the geometry of the 
initial jet profile is evaluated here with a realistic Woods-Saxon distribution, rather than 
using a hard sphere profile. However, these results do show that the photon w 2 is rather 
sensitive to early time dynamics in relativistic nuclear collisions: precise measurements 
thus have the potential to stringently constrain the time-evolution scenarios. 

A global picture is emerging where the physics of QCD jets and the generation of 
electromagnetic radiation in nuclear collisions need to be treated simultaneously and 
consistently. The results displayed here are being generalized to include energy loss by 
elastic collision processes and the effects of a 3D+1 hydrodynamical evolution. 
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